Observational constraints on the dark energy and dark matter mutual coupling 
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Q-c Abstract 

O 

^H ' We examine different phenonienological interaction models for Dark Energy and Dark Matter by perform- 

C^ . ing statistical joint analysis with observational data arising from the 182 Gold type la supernova samples, 

fvq , the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave 

>, . 

^^ I Anisotropy Probe observations, the baryon acoustic oscillation measurement from the Sloan Digital Sky 



o 



X 



Survey and age estimates of 35 galaxies. Including the time-dependent observable, we add sensitivity of 



■^^ ' measurement and give complementary results for the fitting. The compatibility among three different data 

O 

•^ ' sets seem to imply that the coupling between dark energy and dark matter is a small positive value, which 

^ . satisfies the requirement to solve the coincidence problem and the second law of thermodynamics, being 

compatible with previous estimates. 

PACS numbers: 98.80. Cq, 98.80-k 



Our universe is undergoing an accelerated expansion driven by a so called dark energy (DE). Since DE 
occupies almost 70% of the energy content of the universe today, it is natural to consider its interaction 
with the remaining fields of Standard model and its generalizations. It has been claimed that the coupling 
between DE and dark matter (DM) can provide a mechanism to alleviate the coincidence problem [1|, \^ . 
Furthermore, it has been argued that an appropriate interaction between DE and DM can influence the 
perturbation dynamics and affect the lowest multipoles of the CMB spectrum [J, |j] . Recently, it has been 
shown that such an interaction could be inferred from the expansion history of the universe, as manifested 
in the supernova data together with CMB and large-scale structure[5]. However, the observational limits 
on the strength of such an interaction remain weakQ. Signature of the interaction between DE and DM in 
the dynamics of galaxy clusters has also been analyzed [3, Q . Other discussions on the interaction between 
dark sectors can be found in [3, [lO, [lH ■ 

The interaction between DE and DM is a major issue to be confronted in studying the physics of DE. 
However, since neither DE nor DM is actually known to us, it is hard to describe the interaction from 
first principles. Some att emp ts to discriminate the interaction from the thermodynamical point of view 
have been raised recently '12,13]. Most studies on the interaction between dark sectors rely either on the 
assumption of interacting fields from the outset |14l.ll5l|. or from phenomenological requirements [3. 6. 91. Il6l|. 
In view of continuity equations, the interaction between DE and DM must be a function of the energy 
density multiplied by a quantity with units of inverse of time, which can be chosen as the Hubble factor 
H. There is freedom to choose the form of the energy density, which can be any combination of DE and 
DM. Thus, the interaction between DE and DM could be expressed phcnomcnologically in forms such as 

Q = QiHpDAiMM, Q = Q{Hpde)M, Q = Q{H{PDE + Pdm)MM- 

It is of great interest to investigate effects of different forms of interaction between DE and DM on the 
universe evolution. In [l6| the impact of the interaction proportional to the DM energy density on the 
determination of a redshift dependent DE equation of state (EOS) and on the DM density today has been 
studied from SNIa data. It has been shown that the presence of such a coupling increases the tension between 
the CMB data from the analysis of the shift parameter and SNIa data for realistic values of the present DM 
density fraction. Recently, a statistical joint analysis by using observational data coming from the new 182 
Gold SNIa samples, the shift parameter of the CMB given by the three-year WMAP observations and the 
baryonic acoustic oscillation measurement from the SDSS has been carried out for the interaction between 
DE and DM[I3|. Comparisons concerning the influence on cosmological parameters and the effect on solving 
cosmic coincidence problem among different forms of phenomenological interaction models have been done. 
It was argued that consequences of DE and DM interaction on cosmological parameters are sensitive to the 
DE EOS. Choosing an appropriate EOS and the interaction in proportional to the energy density of DE, a 
positive coupling turns out to be more probable and the coincidence problem gets alleviated. However, for 
other forms of the phenomenological interaction models and for other DE EOS, one gets a negative coupling 
between dark sectors which will result in unphysical situations and fail to solve the coincidence problem. The 
negative coupling has also been seen by using the same data from SNIa together with CMB and large-scale 
structure for the interacting holographic DE model with the interaction proportional to the total dark sector 
energy density [5( and other models describing the interaction in proportional to the DM energy density [6|. 
It was argued that the negative coupling is not able to alleviate the coincidence problem[l^ and the model 
does not obey the second law of thermodynamics [ij]. Using the galaxy cluster data, the coupling was 
obtained to be positive indicating the energy decay from DE to DM. 

To reduce the uncertainty and put tighter constraint on the value of the coupling between DE and 
DM, new observables should be added. Recalling that the test of cosmological models by SNIa data is a 
distance based method, it is of interest to look for tests based on time-dependent observables. In [4, 19], the 
age of an old high redshift galaxy has been used to constrain the model. In this work we will combine four 
fundamental observables including the new 182 Gold SNIa samples, the shift parameter of the CMB given by 
the three-year WMAP observations, the baryon acoustic oscillation (BAO) measurement from the SDSS and 
age estimates of 35 galaxies provided in[26| to perform the joint systematic analysis of the coupling between 
dark sectors. We expect that sensitivities of measurements of different observables can give complementary 
results on the coupling between dark sectors. We will compare the compatibility among SNIa data including 
BAO, CMB and age data and determine the tendency of the coupling results. 

Concerning energy conservation for the overall energy density of dark sectors, we can suppose that the 
interaction between DE and DM is described by 

Prn + iHp„i = Q, (1) 

Ph+3H{1+UJD)PD = -Q, (2) 

where Q denotes the interaction term. 



From the equations above, phenomenological forms of the interaction between DE and DM must be a 
function of the energy densities multipHed by a quantity with units of inverse of time, which have possible 
expressions, such as (1) Q = 5H{pdm + Pde), (2) Q = 5HpuM and (3) Q = SHpoE etc. 

We will constrain the coupling between DE and DM in different phenomenological interaction models by 
using the latest observations (golden SN la, the shift parameter of CMB and the BAO) and combining them 
with the lookback time data. We will not specify any special model of DE. Considering recent accurate 
data analysis showing that the time varying DE gives a better fit than a cosmological constant and in 
particular, DE EoS can cross —1 around z = 0.2 from above to below[20j, we will employ two commonly 
used parameterizations in our work, namely 

uJi{z)^wo + jY^^, (3) 

uJii{z)=wo + -^^-^. (4) 

The up-to-date gold SN la sample was compiled by Riess et al [2l|. This sample consists of 182 data, 
which gives the distance modulus at redshift z. The distance modulus is defined as 

/■^ dz' 

Pth{z;P,M) = 5\og,o{dL{z)/Mpc) + 25 = 5\og,o[{l + z) / _-] + 25 - SlogioiJo, (5) 

Jo ^\z ) 

where the luminosity distance dL{z) = ^^ ^ /g -^Itt, the nuisance parameter M = 5 log^g -^^o is marginal- 
ized over by assuming a flat prior P(Hq) = 1 on iJo and P describes a set of parameters characterizing the 
given model. 

An efficient way to reduce the degeneracies of the cosmological parameters is to use the SN la data in 
combination with the BAO measurement from SDSS [22| and the CMB shift parameter [23|]. Using a large 
sample of 46748 luminous red galaxies covering 3816 square degrees out to a redshift of z = 0.47 from the 
SDSS, Eisenstein ct al [2^ have found the model independent BAO measurement which is described by the 
A parameter 

A = ,/n7nEizsAo)-'/n-^ / 7^]^/^ = 0-469(Tf^)-°-^^ ± 0.017, (6) 

zbao Jo ^\z ) U.y» 



where Ug can be taken as 0.95 [2j] and zbao — 0.35. 
The CMB shift parameter is given by 



(7) 



where zis = 1089. This CMB shift parameter R captures how the Z-space positions of the acoustic peaks 
in the angular power spectrum shift. Its value is expected to be the least model independent and can be 
extracted from the CMB data. The WMAP3 data 111 gives R = 1.70 ± 0.03 [11^_^ 

We now turn to the lookback time observations. They have been used in [25| and have been shown 
effective to provide a complementary test of different models. By assuming the total age of the universe to 
be to = 13.7 ± 0.26*2/7% as given by current CMB measurement [2J|, we transform the age estimates of 35 
galaxies provided in [2^. The lookback time-redshift relation is defined by 



dz' 

(l+^0^(z')' 



tUz;P)^H,' I j^-^^^^,^, (8) 



where Hq^ = 9.78h~^Gyr. We have adopted the recent value 0.72 for h given by the HST key project |27| . 
P stands for the model parameters. To use the lookback time and the age of the universe to test a given 
cosmological model, we follow [2^ and consider an object i whose age ti{z) at redshift z is the difference 
between the age of the universe when it was born at redshift zp and the universe age at z, 

Using the lookback time definition, wc have t{zi) — tLizp) — thiz). Thus the lookback time to an object at 
Zi can be expressed as 

tt{z,) = tLizp) - t{z,) = [i°^^ - t,{z)\ - [tf^ - tL{zp)\ = tf-' - t,{z) - d/, (10) 



where df = t°^^ — tL{zF) is the delay factor. 

We employ the Monte-Carlo Markov Chain (MCMC) method [23| to explore the parameter space. By 
using two parameterizations for the EoS of DE iOi,LUji, we show in table 1 the parameter space when the 
coupling between DE and DM is taken proportional to energy densities of DM, DE, total DM plus DE (T), 
respectively. 

Comparing with the result obtained in [17[, it is interesting to find that by adding the new observable, 
the lookback time data, it is possible to have positive coupling between DE and DM, especially for the EOS 
with the form luh. 

Table 1. 

PARAMETERS AT 68.3% CONFIDENCE LEVEL 



Coupling 


EoS 


Wo 


Wl 


^m 


5 


T 


LOI 


J^-^'J-o.is 


1 99+0.20 
-•^•^^-0.27 


r, r,y+0.02 

'-'•^'-0.02 


'J-'J^-0.04 


T 


UJII 


— i.OU_Q 3Q 


3.901^,:°? 


'J-^D-0.02 


Qr,-|+0.03 

^•'^-'^-0.03 


DM 


OJJ 


1 -17+0.16 
^^•J^'-0.14 


-I 9O+0.21 
-'^■^°-0.32 


'^■^'-0.02 


-n n2+°-3i 


DM 


UJII 


-, r;n+0.32 


o 01+212 
"^■^^-2.34 


r, r,r-+0.02 
"■^^-0.02 


O-Olt^ot 


DE 


Wl 


^•^^-0.16 


1 iq+0.19 
-'^•^^-0.28 


n 07+0.01 

'J-^'-o.oi 


-o.o4l°:}| 


DE 


uin 


1 /iQ+0.31 


o ^70+2.13 
'-*• ' °-2.39 


^■^0-0.02 


0.051°:?^ 



Note — (l)The fiat priors on the parameters for EoSl are, -10 < wq < 10, -10 < «)i < 10, < Sim < 0.8, 
-0.5 < 5 < 0.5;(2)For EoS2 ,-10 < too < 10,-15 < wi < 15, < fl,„ < 0.8, -1 < 5 < 1. The CMBEASY GUI is utilized to 
process the MCMC chains. 

Table 2. 
PARAMETERS AT 68.3% CONFIDENCE LEVEL 



Coupling 


EoS 


Wo 


n,n 


s 


T 


Ll>i 


1 10+0.02 
^■^■-'-0.08 


"■^D-0.01 


n2+°-i* 

^■^^-0.03 


T 


LOIJ 


1 cn+0.07 
^•O'J-0.08 


n. r,r-+0.02 

"•^^-0.01 


'-'•'J^-0.02 


DM 


UJI 


1 99+0.04 
^^•^^-0.06 


'J-^O-o.oi 


'-'•"*-0.02 


DM 


LOlJ 


1 fin+007 
^•'3'J-0.08 


"•^^-0.02 


o.oiiro^ 


DE 


UJI 


1 1 D+0.04 
-J^-J^8_o.l2 


0.261°:°} 


0.061°:^;^ 


DE 


LUJI 


1 ^0+0.08 


0.26t°:°i 


0.051°:°^ 



Note — The flat priors are, -5 < u)o < 5, < Sim. < 0.5, -0.2 < <5 < 0.2. 

We perform the data comparisons for different phcnomcnological interaction models between DE and 
DM with two different parameterizations of DE EOS. By taking priority of wi obtained from MCMC fitting 
as central values, we plot the contours in ^^dm — wq plane for different interaction models in Figs 1, 2, 3 
with DE EOS llij. Similar results for DE EOS uijj have also been observed. Green lines indicate the result 
from CMB shift parameters and blue lines are from lookback time result. 

We observe that when |<5| is over the range in table 1, there appear poorer compatibility among the 
three data sets, especially between CMB and SNIa data. For the big positive S, this incompatibility for 
the interaction between dark sectors proportional to the DM energy density was also observed in [iy|. For 
small enough \S\, the range obtained from CMB shift parameters will change more compared to constraints 
from other two data sets. The lower green line for the CMB shift moves upper when small S becomes more 
positive and the upper green line becomes more flattened. This leads more overlapped region for three data 
sets in the floM — wo plane for small positive value of the coupling. However when the positive coupling 
is over a limit, the lower green line will cut the contour, while the upper green line cannot efhciently move 
upper, which will reduce the overlapped region of the constraints from three different data sets. This result 
holds for all forms of phcnomcnological interaction models and different parameterizations of DE EOS. 

Besides, in the small \d\ range, when 6 becomes more positive, we observed that there are more overlaps 
between constraints from the SNIa+BAO and lookback time data sets. Thus from the compatibility of three 
different data sets we obtain the tendency of small positive coupling between DE and DM. 

Choosing now the priority of wi as the central value from MCMC, we obtain the parameter space listed in 
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Figure 1: Contour plots in the Qm-wo plane for a variable EoS uij when Q is proportional to the total energy density of DM 

and DE after giving a prior to wi = 1.22. The compatibility among CMB , SNIa+BAO, Lookback time can be compared by 

examining the 2a contour of CMB shift constraint (green line), the Lookback time constraint (blue line), and Icr, 2ct contours 

of SNIa+BAO rcsult(red line). The 5 years WMAP results for uiq are also indicated by parallel lines. 



table 2. Using the best-fit results of these parameters, we study the coincidence problem. We pay attention 
to the ratio of energy densities between DE and DM, r — pdm / Pde, and its evolution. In Fig. 4 we show 
the behavior of r for the interaction between dark sectors in proportional to DM energy density when we 
choose DE EOS to be w//. For other interaction forms and for DE EOS in the form oiidi^ojn, r behaviors 
are similar. We see that with the positive coupling obtained from the best-fit leads to a slower change of r 
as compared to the noninteracting case. This means that the period when energy densities of DE and DM 
are comparable is longer compared to the noninteracting case. Thus it is not so strange that we now live 
in the coincidence state of the universe. In this sense the coincidence problem is less acute when compared 
with the case without interaction. Similar argument has also been given in [30|. 

It is also worthwhile commenting on the results in the light of the recent 5 years WMAP results |31j . 
We included, in figures 1 and 3 the limits for tuo from 5 years WMAP results. They turn out to be perfectly 
compatible with the SNIa+BAO countours at la. This implies confidence in the results of the present 
paper. For positive coupling, it has more possibility for the overlapped region among three data sets to 
accommodate wq within the 5 years WMAP region, which gives further strength to the claims concerning 
the sign of the interaction. 

In summary, we have examined different phenomenological interaction models between DE and DM by 
performing statistical joint analysis with observational data from the new 182 Gold SNIa samples, the shift 
parameter of the CMB given by the three-year WMAP, the baryon acoustic oscillation measurement from 
the SDSS and age estimates of 35 galaxies. Comparing with the test by just using data from SNIa together 
with CMB and large-scale structure [13], we observed that adding the age constraint, we get a tendency 
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Figure 2: Contour plots in the Qm-WQ plane for a variable EoS a;/ when Q is proportional the energy density of DM after 

giving a prior to wi = 1.28. Green lines are for 2a contours of CMB shift constraint, blue lines are for the Lookback time 

constraints and red lines are for la ,2a contours of SNIa+BAO results. 



towards a positive coupling between DE and DM, especially for the DE EOS with the form oju. This shows 
that the new observable can add sensitivity of measurement and give a complementary result for the fitting. 
We have studied the compatibility among three different data sets including SNIa plus BAO, CMB 
shift and lookback time. We found that the bigger couplings \6\ between dark sectors lead to a poorer 
compatibility, especially comparing CMB with other two data sets. For small \5\, we observed, for all 
phenomenological forms of interaction with two parameterizations of DE EOS, the same tendencies for the S 
to be a small positive number. The small positive coupling result is consistent with that got independently by 
galaxy cluster analysis [8| . The positive coupling is required to alleviate the coincidence problem and avoid 
some unphysical problems met in [la . |l7| . It is also the requirement of the second law of thermodynamics 
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Figure 3: Contour plots in the Om-'i'o plane for a variable EoS uij when Q is in proportional to the DE energy density after 

giving a prior to wi = 1.19. Green lines are for 2a contours of CMB shift constraint, blue lines are for the Lookback time 
constraints and red lines are for Icr ,2a contours of SNIa+BAO results. The 5 years WMAP results for wq are also indicated 

by parallel lines. 




Figure 4: The red line indicates the evolution of the ratio of energy densities between DE and DM when the interaction is in 
proportional to the energy density of DM and DE EoS is in the form of cou. We have compared the interacting case(red line) 

with the non-interacting case(blue line). 
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